Correlation functions and queuing phenomena in growth processes with drift 
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We suggest a novel stochastic discrete growth model which describes the drifted Edward- Wilkinson 
(EW) equation dh/dt = ud 2 h — vd x h+ri(x, t). From the stochastic model, the anomalous behavior of 
the drifted EW equation with a defect is analyzed. To physically understand the anomalous behavior 
the height-height correlation functions C(r) = {\h(xo + r) — h(xo)\) and G(r) = (\h(xo + r) — h(xo)\ 2 ) 
are also investigated, where the defect is located at xo. The height-height correlation functions follow 
the power law C(r) ~ r° and G(r) ~ r a with a' = a" = 1/4 around a perfect defect at which no 
growth process is allowed, a — a" = 1 /4 is the same as the anomalous roughness exponent a = 1 /4. 
For the weak defect at which the growth process is partially allowed, the normal EW behavior is 
recovered. We also suggest a new type queuing process based on the asymmetry C(r) ^ C(—r) of 
the correlation function around the perfect defect. 

PACS numbers: 68.35.J, 05.10.G, 11. 10. J 



I. INTRODUCTION 

Occasionally some of the research results on older top- 
ics give excellent insights to the subsequent researches. 
One of the such works is the recent work on the drifted 
Edward- Wilkinson (EW) equation 



dh(x, t) 
dt 



vd%h(x, t) — vd x h(x, t) + rj{x, t) , (1) 



with a fixed boundary condition (FBC) h(xo — 0, t) — 
h(xo = L,t) = pj. If the diffusion term ud^.h(x,t) is 
ignored in Eq. Q, one can obtain roughness exponent 
a = 0, growth exponent (3 — and dynamic exponent 
z = 1 by a scaling argument and some other methods 

SB- But for periodic boundary condition (PBC), 
the drift term —vd x h(x, t) is irrelevant, and thus Eq. (JTJ 
becomes physically equivalent to the normal EW equa- 
tion Hi 



dh(x, t) 
dt 



vd 2 x h(x,t) + ri(x,t) 



(2) 



with a = 1/2, = 1/4, and z = 2. In contrast, if FBC in- 
stead PBC is imposed, the roughness of surface described 
by Eq. Q is exactly proved to show an anomalous be- 
havior with the exponents 0, 



= 1/4, 0=1/4, z = l 



(3) 



The exact solution for Eq. (|TJ 1] gives us an excellent 
insight to understand the role of the drift term and the 
boundary condition as a relevant extension of the normal 
EW equation. 

The drift term and the diffusion term are crucial for 
anomalous exponents under FBC p|. Without the drift 
term, the maximum time which the roughness is fully 
developed for the system is order of L 2 jv. If the sys- 
tem is affected by the drift, the time scale that any 



noise-generated local structure starts from one bound- 
ary to the other boundary is t x = L/v. The fluctu- 
ation of the surface width is saturated after this time 



scale. Therefore W 2 



1/2 
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L 1 ! 2 and thus one can 
get the anomalous exponents (J3J. In Ref. l], local 
width is numerically studied to test the validity of this 
physical mechanism. The local width is increased as 
W 2 (xi) oc i; 1 / 2 ~ {xi/v) 1 / 2 , with xi being the position 
where the width is measured. But the definition of the 
local width is not physically clear and there is a sort of 
difficulty to measure W 2 (xi) for the small system. 

For last two decades, the stochastic discrete growth 
models whose dynamical scaling behavior is identical to 
those of a given continuum growth equation have been 
used to investigate the dynamical scaling behavior for 
the corresponding surface roughness 0, U IE E 1 ] ■ The 
stochastic discrete growth model is a powerful tool which 
simplifies the complicate growth behavior and provides 
an essential link between theory and experiments. For 
example, the surface growth which is described by nor- 
mal EW equation has been studied through the Family 
model 10]. But no simple growth model which exactly 
corresponds to the drifted EW equation with FBC has 
been suggested yet. Other existing models cannot ex- 
plain the anomalous behavior of the surface fluctuation 
which comes from the drift term and FBC in Eq. Q). 
First motivation of this paper is thus to make a simple 
stochastic discrete model which satisfies Eq. Q with 
FBC using a stochastic analysis method 0] ■ 

Second motivation of this paper is in the explanation 
of the physical mechanism which causes the anomalous 
behavior by using a height-height correlation func- 
tion. The height-height correlation function can be easily 
treated through our handy stochastic discrete model. As 
mentioned above, the local width Q] is somewhat imma- 
ture physical quantity and has not been widely used for 
the analysis of surface roughness 0. In contrast, height- 
height correlation function has been extensively used to 
show the exact dynamic scaling in the surface growth 
physics 6] . The height-height correlation function which 
we will use to analyze the surface with FBC is an aver- 
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age of height differences or of squared height differences 
between h(xo + r,t) and h(xQ,t), where r (1 < r < L) 
is the distance from a boundary xq — at time t. If we 
know the height-height correlation function between two 
fixed boundaries, we can get clean and clear information 
how the saturated surface of whole system is formed. We 
will measure the height-height correlation function and 
explain the anomalous behavior J3j through it. 



Third motivation of this paper is to suggest about a 
new type queuing phenomena originated from the anal- 
ysis of the height-height correlation function. Queuing 
|l2j is a common nonequilibrium phenomenon in nature. 
It is well established that many driven flow processes 
belong to the same universality class as Kardar-Parisi- 
Zh ang (KPZ) type growth of one dimensional interfaces 
0, El- For example, a traffic jam [l4| caused by slow 
and fast bond of the so-called asymmetric simple exclu- 
sion process (ASEP) ^| is related to the faceting on 
KPZ growth. ASEP breaks translational invariance in 
many ways. One of them is to retain periodic boundary 
conditions, but to introduce a defect into the system by 
modifying the transition rates locally [T3 . ITfil | . By prop- 
erly chosen injection/removal rates or defect strength, 
nonequilibrium phase transitions occur between profiles 
of different shapes and average densities. In our model, 
we map the fixed boundary condition to a defect which 
exists on the center of system. The defect site means that 
the growth at a certain site is impossible or is controlled 
by a probability p which is called a defect strength. That 
is, the fixed boundary condition h(xo = 0, t) — h(xo = 
L, t) — is shifted by the mapping xq — > x + L/2 to the 
defect site h(xo = L/2,t) = with the periodic bound- 
ary condition. For p — any deposition or evapora- 
tion process at xq cannot occur (h(xo) = 0), since the 
strength of the defect is absolutely strong and thus the 
surface heights cannot be freely increased by the defect. 
If p =/= 0, the deposition or evaporation process at x$ is 
accepted by p and the system with p = I is the same as 
PBC without the defect. The height-height correlation 
function of our stochastic discrete model shows lateral 
asymmetry C(r) C(— r) around the defect, in contrast 
to KPZ-typc queuing phenomena or ASEP which shows 
up-down symmetry breaking under the transformation 
h(x) — > —h(x). The new queuing phenomena from the 
drifted EW equation with a defect is physically dif- 
ferent from usual ASEP or queuing phenomena which 
is related to the nonlinearity of KPZ equation 0, be- 
cause it is from the linearity of the drifted EW equation 
l(T)l. Therefore it is very interesting to compare our new 
queuing phenomena to earlier studied KPZ type queuing 
phenomena. In this paper, we want to investigate new 
queuing phenomena from the defect and the drift using 
our stochastic discrete model and apply it to the parking 
garage model [l4| . 



II. A NOVEL METHOD TO MAKE 
STOCHASTIC DISCRETE MODEL 

In this section, we introduce a new method to make a 
discrete stochastic growth model which physically corre- 
sponds to the continuum equation. Normally the surface 
configuration is described by the continuum Langevin 
equation |6( 



dh(x, t) 

at 



J 7 (h(x),t)+r)(x,t) , 



(4) 



where !F{h{x),t) depends on the surface configuration 
{h(x)}. In the method, the surface configurations are 
described in terms of integer height variables as H = 
{hijfLi — 0' ^2, • ■ • on a lattice. Then the continuum 
equation can be discretized into 



ot 



(5) 



Associated with Eq. JSJ, the general master equation 
|l7l Il8l Il9j ] for the surface evolution can be written as 



>(H\ H)P(H', t)-J2 H')P(H, t) 



dP(H,t) 

Ft 

//' It' 

(6) 

Here uj{H' ,H) is the transition rate from the configura- 
tion H to the configuration H' and P(H, t) is the prob- 
ability that the system is in the configuration H at time 
t. Using the transition moments |19| 



Kf\H) = Y\{h' i -h i )uj{H',H) 



H 



and 



K?\H) = 5>J - hj)w(H',H) 



(7) 



(8) 



H' 



the master equation can be rewritten as Fokker- 
Planck equation [l7LliaTl9| 



dP(H,t) 
dt 



1 d 2 



^ 2 dhidhi 13 

h3 



[K%\H)P(H,t)], (9) 



by Kramers-Moyal expansion. Therefore we can obtain 
the discrete version of equivalent Langevin equation from 
Eq. © as 



(10) 



and thus is equal to Ti of Eq. JSJ. Here the 

noise Tji(t) has the properties < rji(t) >— and < 
mitfrjjit') >= K^S^ - t'). When the noise is white, 

(2) 

K\- ' is equal to 2DSij . 
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Normally the evolution process in a stochastic surface 
growth model is defined by one-particle deposition or 
evaporation in a unit evolution process 6] . Hence if only 
one-particle deposition (hi — > hi+a) or evaporation^ — » 
hi — a) is imposed at a randomly chosen site i in a unit 
evolution process, then k\ 1] and are 



= v 2 [h 



i+l 



hi-i - 2hi] - -[h i+1 



I -(14) 



Kf'(H) = auj id - aui, 
K$\B) - 







a 2 u> id + a 2 uj ie ,i = j 



(11) 



from Eqs. J7J) and ©, where uJid is 

u>({hi, h 2 , ■ ■ ■ ,hi,--- , h N }, {hi,h 2 , ■■ ■ , h + a, ■ ■ ■ , h N }) 
and uj ie is u({h 1 ,h 2 , • • • , hi, ■ ■ ■ , h N }, {hi,h 2 , ■■■ ,h - 
a, ■ ■ ■ , /iAr}).Then from the Eq. these rates are 



"«= + •© + - (i 2 ) 




(13) 



Here we impose the unit step function &(x), i.e. O(ir) = 1 
for x > and O(x) = for x < 0, to the transition rates 
tOid and LUie, because the transition rates must be always 
positive. If r]i(t) is the local white noise, the evolution 
site i is randomly chosen to apply the rates (12) and (13). 

Now we want to describe the discrete stochastic 
model which exactly correspond to Eq. Q on an one- 
dimensional (Id) substrate in detail. Generalization to 
that on higher-dimensional substrates can easily be ob- 
tained from the description of the Id model. The discrete 
model is defined by the following steps, (i) Select a site 
i randomly, (ii) Calculate the deposition rate (the rate 
for hi — > hi + 1) uiid and the evaporation rate (the rate 
for hi — > hi - 1) u ie from Eqs. ifl^l and JI3J). (We set 
o = l.) (iii) Compare a random number R (0 < R < 1) 
with the deposition probability Pid = u>id/(ujid + ^>ie) 
and the evaporation probability Pi e = uji e / (uji d + u>i e )- If 
R < Pid, take the deposition process. Otherwise take the 
evaporation process. 



III. DYNAMICAL SCALING PROPERTIES OF 
THE STOCHASTIC MODEL 

We now report the numerical results of our stochas- 
tic discrete model for the drifted EW equation with the 
defect which is controlled by a strength p. The initial sur- 
face is always set to be fiat (hi(t — 0) = 0) and the lateral 
periodic condition is always imposed. All data are taken 
by average over more than 500 independent runs. We 
first explain the scaling properties of the surface width 
for several defect strengths p. The first transition mo- 
ment K^- 1 ' of drifted EW equation (JIJ at a randomly 
selected site % is 



We use the parameters v% = 1, v = 1 and D = 0.01 for 
Eq. (JTJ. We check the dynamical scaling behaviors for 
other parameter values to obtain the same behaviors for 
u 2 = 1, v = 1 and D = 0.01. Used system sizes are 
L = (2 6 + 1), • • • , (2 12 + 1) and the defect is located at 
the center site xo = (L — l)/2 + 1 of all systems. 

Usually the surface roughness is defined as W(L,t) = 



(h 2 ) — (h). The scaling relation of W satisfies the 
following behavior [(| 



W(L,t) 



L a f(t/L z ) 

, t < L z 
L a , t > L z 



(15) 



K, 



(i) 



[u 2 V 2 h - vVh]i 



Figure 1(a) shows the obtained saturated surface width 
W sa t(— W(t 3> L z )) from the stochastic model for vari- 
ous system size L and defect strength p. W sa t satisfies the 
scaling relation W sat ~ L a with a = 0.50(1) for p = 0.1 
and p = 1. Since p = 1 case is the same as that with- 
out the defect, a should be EW value a = 1/2. Even 
though the defect strength is sufficiently large (p = 0.1), 
the surface roughness still follows normal EW behavior. 
In contrast, for perfect defect or p = 0, the surface width 
W follows the anomalous behavior with a — 0.25(1) in- 
stead of EW behavior. The main plot of Fig. 1(b) shows 
that W(L,t) for various L and t collapses well to scaling 
function Eq. 1|15|) for p = 0.1. We use the exponents 
a = 0.50 from the fitting of W - L a (t > L z ) in Fig. 
1(a) and j3 = 0.25 (W ~ t^(t < L z )) as shown in the 
inset of Fig. 1(a). Therefore the dynamic exponent is 
z = a/P = 2 which is the value of the normal EW uni- 
versality class. The results in Fig. 1(b) show that the 
dynamical scaling property for p ^ satisfies the normal 
EW behavior. In contrast, the exponents are changed for 
p = (Fig. 1(c)). The time dependence of W(L,t) of 
the model at p = is shown in the inset of Fig. 1(c). 
[3 = 0.249(1) in Fig. 1(c) is obtained by applying the 
relation W ~ t (t < L z ) to W(t) of the system with 
L = 2 12 + 1. We also check the dynamical scaling rela- 
tion of Eq. by plotting ln(W/L a ) against ln(t/L z ) 
in Fig. 1(c). The data collapse well to the scaling func- 
tion Eq. I|15l) with a — 0.25 and 2 = 1. The results in 
Fig. 1 clearly show that the anomalous behavior occurs 
only when the perfect defect (p = 0) exists in the system. 
These results agree with the analytic result of Ref. pj. 
From these results, we confirm that a sudden crossover 
from the anomalous behavior to the normal EW behavior 
occurs as soon as p deviates from p = 0. 

To understand the physical origin of the anomalous 
interface profile which is induced from the defect (FBC) 
and the drift term in Eq. (Q, we study several height- 
height correlation functions. The height-height correla- 
tion function which has been extensively used for the 
analysis of surface roughness 6] should be one of the effi- 
cient methods to see the physical mechanism. Normally 
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FIG. 1: (Color online) (a) Plot of InVKsat against InL in 
the saturation regime (t 2> L z ). The straight lines rep- 
resent the relation W sa t ~ L a . a = 0.25(1), 0.5(1), 0.5(1) 
for p = 0, p — 0.1, and p = 1, respectively. Used system 
sizes are L = (2 6 + 1), ■ ■ ■ , (2 12 + 1) and used parameters are 
v = 1, v = 1, £> = 0.01. (b) The plot of ln(W/L a ) for the dif- 
ferent L against ln(t/L z ) with a — 1/2 and z — 2 at p — 0.1. 
The numerical data for VK nicely collapse to one curve which 
satisfies the relation 1151 . (c) Scaling plot showing that the 
data for \n(W/L a ) plotted against \n(t/L z ) for various L col- 
lapse to a single curve supporting the scaling function 115H 
with a = 1/4 and z = 1 at p = 0. Insets of (b) and (c) are 
log-log plots of W against t. 
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FIG. 2: (Color online) (a) Plots of the height-height corre- 
lation function 117II for p = and p — 0.1 when t ^> L z . 
Inset shows the log-log plot of the dotted square part in 
the main plot. C(r) satisfies the relation C ~ r a with 
a' = 0.55(1) (p = 0.1) and a' = 0.25(1) (p = 0). (b) Plots of 
the height-height correlation function Eq. 1181 for p = and 
p = 0.1. Inset shows the log-log plot of the dotted square part 
in the main plot. G(r) satisfies the relation G ~ r 2a with 
2a" = 1.10(1) (p = 0.1) and 2a" = 0.50(1) (p = 0). Used 
system size is L = 1025. 



the height-height correlation function 6] means 



G 



PBC 



i L 

-^([h(x,t)-h(x + r,t)] 2 ) , 



(16) 



x = l 



when the periodic boundary condition is imposed. How- 
ever when a local defect is the relevant factor for the 
dynamical scaling properties of the interface, the cor- 
relation function which measures the height differences 
between h{xo) and h(x,t)(x ^ xq) should be more pow- 
erful one. One of such correlation functions is the naive 
and simple height-height correlation function C (r, t) = 
(h(xo +T,t) — h(xQ, £)). Here ( ) means the average over 
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independent runs and r means the position from the de- 
fect Xo(= L/2). However, C(r, t) cannot characterize the 
surface configuration physically well for the equilibrium 
surfaces in which the deposition and evaporation pro- 
cesses occur with nearly the same probability. The equi- 
librium surface fluctuates around (h) = and C(r, t) has 
no information for the surface height correlation around 
the defect. Instead, we use two height-height correlation 
functions C(r, t) and G(r, t) defined as 



(a) 



and 



C(r,t) = (\h(x +r,t)-h(x ,t)\) (17) 



G(r,t) = (\h(x +r,t)~h(x ,t)\ 2 ) , (18) 



respectively. For the analysis of the saturated surface 
configurations, we concentrate only on the saturated cor- 
relation functions C(r,t ^> L z )(= C(r)) and G(r,t ^ 
L z )(= G(r)). Used system size for the measurement of 
C(r) and G(r) is L = 1025 and thus the defect site is 
located at xo = 513. These two height-height correla- 
tion functions are expected to scale as C(r) ~ r a and 



G(r) 



J2a" 



as shown in Figs. 2(a) and (b). Furthermore 



for p — we obtain the exponents a 1 = a" — 0.25(1), 
which is very close to the roughness exponent a — 0.25(1) 
in Fig. 1(a). For p = 0, 



and 



C(r) = (\h(x + r,t)\) 



G(r) = (\h(x +r,t)\ 2 ) 



(19) 



(20) 



and thus we clearly see the (\h(xo + r,t)\) ~ r 1 / 4 , which 
decide the anomalous behavior a — 1/4. This means 
the saturated surface configurations for the drifted EW 
equation with a defect satisfies 



W 



(\h(x +r,t)\) 



x 



r,t)) 



(21) 



Therefore physically the surface configuration is decided 
by the FBC (or the defect), which does not allow the sur- 
face fluctuation freely. However, PBC allows the surface 
fluctuation freely even though there exists the drift term 
d x h. 

In contrast, as soon as the growth at the defect is 
allowed even with a very low probability, we get a' = 
a" w 0.5, which is almost the same as roughness expo- 
nent a = 0.5 of EW universality class. (See the data for 
p = 0.1 in Fig. 2.) Therefore the anomalous behavior 
a = a' = a" = 1/4 for p — disappears and are changed 
into a = a' = a" = 1/2 as soon as the defect condition 
is released or for p ^ 0. 

Even though the defect exists, the drift term does not 
change the universality of the dynamical surface scal- 
ing for < p < 1. Therefore the height-height corre- 
lation function has a lateral symmetry C(r) = C(— r) or 
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FIG. 3: (Color online) Schematic representation of the queu- 
ing phenomena which emerges in the stochastic growth model 
with a defect. Periodic boundary condition is used, (a) The 
initial distribution of cars, (b) A steady state configuration. 
Arrows denote the local velocities of cars. Cars move out from 
garage with a very low speed and they are packed at the exit, 
because of the defect. In the middle of the road, cars main- 
tain their average velocity. On the contrary, the speed with 
which cars get in the garage is much faster than the average 
velocity. 
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FIG. 4: (Color online) Plot of \d x h\ against the position i. 
Used system size is L — 1025. 



G(r) — G{—r) for p ^ as for the normal EW equation 
(J2J. For p = 0, the noise-generated surface structure from 
the defect is pushed by the drift and transmitted all over 
the system. Then the flow of the surface structure comes 
back through the boundary suddenly disappears at the 
defect. Therefore the height-height correlation function 
shows the asymmetry C(r) ^ C{—r) and G{r) ^ G(—r). 
This asymmetric behavior of C(r) and G(r) is a unique 
behavior for the drifted EW equation with a perfect de- 
fect. In contrast the height-height correlation function 
around a defect for KPZ-type surface growth Q al- 
ways shows the lateral symmetry C(r) = C(—r). In 
KPZ-type surface growth model with the defect [l4[, the 
height- height correlation function C(r,p) changes from 
C(r,p) - A\r\ + B\r\ a ' for p < p c to Cir,p) - B\r\ a ' 
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for p > p c - It means that the queueing effect of the de- 
fect spreads out entire system (jp > p c ) or affects only 
around region of the defect (p < p c ). The phase transi- 
tion argument of KPZ-type surface growth is analogous 
to the queuing effect by a defect which is located at a 
certain point of lane [lj]- On the other hand, the height- 
height correlation functions which we measured always 
maintains the simple power behavior C(r) ~ r a and 
G(r) ~ r a . Instead, it shows the sudden crossover of 
the exponents a' and a" as soon as the defect is weakly 
imposed (or p / 0). That is, the effect of the defect 
and the drift influences only when the defect strength is 
absolutely strong (or p = 0). 

IV. APPLICATION TO QUEUING 
PHENOMENA 

The asymmetric form of the interface profile by the 
defect and drift can be applied to a kind of the queu- 
ing phenomena or car condensation in a parking garage 
model [lj, Ujj. Let's think of the traffic flow of cars 
on an one-dimensional road with a periodic boundary 
condition (or a ring- type road). The defect site in our 
model is then reinterpreted as the parking garage (See 
Fig. 3). Initially all cars on the road are equally dis- 
tributed everywhere on the road as shown in Fig. 3(a). 
If the drift EW equation (JIJ with a defect (garage) de- 
scribes a queuing phenomena, the absolute gradient of 
the height \d x h(t)\i = \hi(t) — /ij_i(£)| at i-th site can 
be mapped to the effective speed Vi of a car around the 
site i at the given moment t due to the drift term 
In the steady state the car density pi at the site i is 
then inversely proportional to \d x h\i = — as 
v t (x \d x h\i oc 1/p. 

Figure 4 shows l^a^li for p = in the saturated state 
(or steady state) when t ^> L z . \d x h\i increases very 
rapidly as i increases from the defect position and im- 
mediately reaches a certain average value. Average value 
remains until i approaches the defect site from the left, 
and | d x h \ i grows very rapidly to the highest value in the 
limit i — > x — . These results correspond to the follow- 



ing steady-state traffic flow. The cars slowly move out 
from the parking garage by a sort obstruction at the exit 
or by a lazy parking charge collector. Thus they are 
packed near the exit of the garage and the density p of 
cars at the exit is very high. As soon as cars enter the 
road, the speed and car density reach to the average value 
|<9 x /i|i = v a ~ 1/ p a ~ const, and maintain the average 
values in the middle of road. When cars come close to 
the entrance of the garage, the speed of cars gets much 
faster. Therefore the density of cars at the entrance be- 
comes very low. 

V. SUMMARY 

In this paper, we presented a simple stochastic dis- 
crete model which describes the drifted EW equation JTJ 
with a defect. From the simple model, we can easily 
show the anomalous behavior of the surface fluctuation 
by the drift and defect. The scaling exponents show the 
anomalous behavior a = 1/4, (3 = 1/4, and 2 = 1, only 
when the defect strength is p = 0. The exponents are 
suddenly changed to normal EW exponents as soon as 
the defect strength is weakly imposed (p ^ 0). We also 
measure the height-height correlation functions to char- 
acterize the asymmetric interface profile and to explain 
the exact physical role of the drift term and defect. The 
height-height correlation functions show that the drift 
and the perfect defect (p — 0) makes the asymmetry 
C( r ) ^ C(-r) or G(r) ^ G{-r). These height-height 
correlation functions satisfy the power law C(r) ~ r a 
and G{r) ~ r 2a " with a' = a" = 1/4. The height-height 
correlation function exponents and the roughness expo- 
nent have the same value a — a' — a" — 1/4. Therefore 
the saturated surface configuration for drifted EW equa- 
tion with the defect is solely decided by the defect. The 
asymmetry of the height-height correlation functions is 
very unique when compared to the correlation functions 
for KPZ-type growth model with the defect which always 
shows the symmetry C(r) — C(—r). From this asymme- 
try we suggest a new queuing process around the defect. 
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